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Abstract 

A mixed-symmetry nuclear shell-model scheme for carrying out 
calculations in regimes where there is a competition between two or 
more modes is proposed. A one-dimensional toy model is used to 
demonstrate the concept. The theory is then applied to '^'^Mg and 
'^^Ti. For lower p/-shell nuclei such as '^^~'^^Ti and '^^Cr there is 
strong SU{3) symmetry breaking due to the spin-orbit interaction. 
However, the quadrupole collectivity as measured by B{E2) transition 
strengths in the yrast band remain high even though SU (3) appears 
to be broken. Some results for the so-called A (5) symmetry that falls 
along the C/(5) <-> SU{'i) leg of the Interacting Boson Model are also 
considered. The results show that the mixed-symmetry concept is 
effective, even when strong symmetry breaking occurs. 



1 The Mixed-Mode Concept 

Two dominant and often competing modes characterize the structure of 
atomic nuclei. One is the single-particle structure that is demonstrated by the 
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validity of the mean-field concept; the other is the many-particle collective 
behavior that is manifested through nuclear deformation. The spherical shell 
model is the theory of choice whenever single-particle behavior dominates [1]. 
When deformation dominates, the Elliott SU{3) model or its pseudo-SU(3) 
extension is the natural choice [2]. This duality manifests itself in two domi- 
nant components in the nuclear Hamiltonian: respectively, the single-particle 
term, Hq = J^i^i^i, and a collective quadrupole-quadrupole interaction, 
Hqq — Q ■ Q- It follows that a simplified Hamiltonian H — ^i^i — xQ " Q 
has two solvable limits. 




Figure 1: A schematic representation of the potential for a particle in one- 
dimensional (ID) box at high energies and a harmonic oscillator (HO) restor- 
ing force for low energies. 

To probe the nature of such a system, one can consider a simpler prob- 
lem: the one-dimensional harmonic oscillator in a box of size 2L [3] . As for 
real nuclei, this system has a finite volume and a restoring force of a har- 
monic oscillator type, For this model, shown in Figure 1, there is a 
well-defined energy scale that measures the strength of the potential at the 
boundary of the box, Ec = a;^L^/2. The value of Ec determines the nature 
of the low-energy excitations of the system. Specifically, depending on the 
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value of Ec there are three spectral types: 

(1) For — s> the spectrum is simply that of a particle in a box; 

(2) At some value of u, the spectrum begins with Ec followed by the spec- 
trum of a particle in a box perturbed by the harmonic oscillator poten- 
tial; 

(3) For sufficiently large u the spectrum is that of a harmonic oscillator 
below Ec which is followed by the perturbed spectrum of a particle in 
a box. 



1.0 




State Number (n) 

Figure 2: The relative deviations from the exact energy eigenvalues for the 
ID box plus HO potential (Figure 1) with uj = 16, L = 7c/2, h/2TT = m = 1. 
The open circles represent deviation of the exact energy eigenvalue from 
the corresponding harmonic-oscillator eigenvalue (1 — Eho/ Eexact), the solid 
diamonds are the corresponding relative deviation from the energy spectrum 
of a particle in a ID box, and the solid squares are the first-order perturbation 
theory estimates using particle in a ID box wavefunctions. 

The last scenario (3) is the most interesting since it provides an example 
of a two-mode system. For this case the use of two sets of basis vectors, 
one representing each of the two limits, has physical appeal, especially at 
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energies near E,.. One basis set consists of the harmonic oscillator states; 
the other set consists of basis states of a particle in a ID box. We call this 
combination a mixed-mode / oblique-basis approach. In general, the oblique- 
basis vectors form a nonorthogonal and overcomplete set. Even though a 
mixed spectrum is expected around Ef.^ our numerical study that includes 
up to 50 harmonic oscillator states below E^., shows that the first order 
perturbation theory in energy using particle in a ID box wave functions as 
the zero order approximation to the exact functions works quite well after 
the breakdown of the harmonic oscillator like spectrum. This observation is 
demonstrated in Figure 2 which shows the relative deviations from the exact 
energy spectrum for a particle in a ID box. 

Although the spectrum seems to be well described using first order per- 
turbation theory based on a particle in a ID box wave functions, the exact 
wave functions near E^. have an interesting structure. For example, the zero 
order approximation to the wave function used to calculate the energy may 
not be present at all in the structure of the exact wave function as shown 
in Figure 3. Another feature also seen in Figure 3 is the common shape of 
the distribution of the non-zero components along the particle in a ID box 
basis. The graph in Figure 4 shows this same effect in nuclei, which is usually 
attributed to coherent mixing [4, 5]. 

2 Fermion-based Applications of the Theory 

An application of the theory to '^'^Mg [6], using the realistic two-body interac- 
tion of Wildenthal [7], demonstrates the validity of the mixed-mode concept. 
In this case the oblique-basis consists of the traditional spherical states, that 
yields a diagonal representation of the single-particle interaction, together 
with collective S'f/(3) configurations, that yields a diagonal representation 
of the quadrupole-quadrupole interaction. The results shown in Figures 5 
and 6 illustrate typical outcomes. For example, a SM(2) + (8,4)&(9,2) model 
space (third bar in Figure 6) reproduces the binding energy (within 2% of 
the full-space result) as well as the low-energy spectrum. For this case the 
calculated eigenstates have greater than 90% overlap with the full-space re- 
sults. In contrast, for a pure m-scheme spherical shell-model calculation one 
needs about 60% of the full space, SM(4) - the fourth bar in Figure 6, to 
obtain comparable results. 
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Figure 3: Bar charts that illustrate the similar, coherent structure of the 
25th, 27th and 29th exact eigenvectors in the basis of a free particle in a ID 
box. Parameters of the Hamiltonian are u = 16, L = 7r/2, /i/27r = m = 1. 
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Figure 4: Bar charts that illustrate the very similar, coherent structure of 
the first three yrast states in ^^Cr calculated using the realistic Kuo-Brown-3 
interaction {KB3). The horizontal axis is C2 of SU{3) while the height of 
each bar gives the contribution of that configuration to the corresponding 
yrast state. 
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Figure 5: The graph shows the calculated ground-state energy for Mg as a 
function of various model spaces. SM(n) denotes a spherical shell model cal- 
culation with up to n particles outside of the 0^5/2 sub-shell. Note the dramatic 
increase in binding (3.3 MeV) in going from SM(2) to SM(2)+(8,4)&(9,2) (a 
0.5% increase in the dimensionality of the model space). Enlarging the space 
from SM(2) to SM(4) (a 54% increase in the dimensionality of the model 
space) adds 4.2 MeV to the binding energy. 
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Figure 6: Bar chart that shows representative overlaps of pure SM(n), pure 
SU(3), and obhque-basis results with the exact full sd-shell eigenstates. A 
number within a bar denotes the state with the overlap shown by the bar if 
it is different from the number for the exact full-space calculation shown on 
the abscissa. For example, for SM(2) the third eigenvector has the largest 
overlap with the fourth exact eigenstate, not the third, while the fifth SM(2) 
eigenvector has greatest overlap with the third exact eigenstate. 
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Studies of the lower p/-shell nuclei ^4-482^^ ^^^^ ^^Cr [4], using the realistic 
Kuo-Brown-3 (KB3) interaction [8], show strong SU{2)) symmetry breaking 
due mainly to the single-particle spin-orbit sphtting. Thus the KB3 Hamilto- 
nian could also be considered a two-mode system. This is further supported 
by the behavior of the yrast band B(E2) values that seems to be insensi- 
tive to fragmentation of the SU{'i) symmetry. Specifically, the quadrupole 
collectivity as measured by the B(E2) strengths remains high even though 
the SU (3) symmetry is rather badly broken. This has been attributed to a 
quasi-SU(3) symmetry [5] where the observables behave like a pure SU{?>) 
symmetry while the true eigenvectors exhibit a strong coherent structure 
with respect to each of the two bases. This provides justification for further 
study of the implications of mixed-mode shell-model studies. 

Next we discuss obhque-basis type calculations for ^^Ti using the KB3 
interaction [8]. We confirm that the spherical shell model (SSM) provides 
a significant part of the low-energy wave functions within a relatively small 
number of SSMC while a pure SU (3) shell-model with only few SU (3) irreps 
is unsatisfactory. This is the opposite of the situation in the lower sii-shell. 
Since the SSM yields relatively good results for SM(2), combining the two 
basis sets yields even better results with only a very small increase in the 
overall size of the model space. In particular, results in a SM(2)+SU(3) model 
space (47.7% + 2.1% of the full p/-shell space) are comparable with SM(3) 
results (84%). Therefore, as for the sd-shell, combining a few SU(di) irreps 
with SM(2) configurations yields excellent results, such as correct spectral 
structure, lower ground-state energy, and improved structure of the wave 
functions. However, in the lower srf-shell S'f/(3) is dominant and SSM is 
recessive (but important) and in the lower p/-shcll one finds the opposite, 
that is, SSM is dominant and 5^7(3) is recessive (but important). 

For ^^Ti adding the SU{3) to the SM(2) increases the model space from 
47.7% to 49.8% and gives results that are slightly better than the SM(3) 
which is 84% of the full space. In '^^Mg the position of the K=2 band head 
is correct for the SU(3)-type calculations but not for the low-dimensional 
SM(n) calculations [6]. In "^^Ti it is the opposite, that is, the SM(n)-type 
calculations reproduce the position of the K—2 band head while SU(3)-type 
calculations cannot. Furthermore, the low-energy levels for the SU{3) case 
are higher than for the SM(n) case. Nonetheless, the spectral structure in 
the oblique-basis calculation is good and the SM(2)-|-(12,0)&(10,1) spectrum 
(~50% of the full space) is comparable with the SM(3) result (84%). 
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The overlaps of SU(3)-type calculated eigenstates with the exact (full 
shell-model) results are, often less than 40%, but the SM(n) overlaps are 
considerably bigger with SM(2)-type calculations yielding an 80% overlap 
with the exact states while the results for SM(3) show overlaps greater than 
97%, which is consistent with the fact that SM(3) covers 84% of the full space. 
On the other hand, SM(2) + (12,0)&(10,l)-typc calculations yield results that 
are as good as those for SM(3) in only about 50% of the full-space and 
SM(1)+(12,0)&(10,1) overlaps are often bigger than the SM(2) overlaps. The 
oblique-bases SM(2) + (12,0)&(10,1) results for *^Ti (^50%) yields results 
that are comparable with the SM(3) results (~84%). In short, the oblique- 
basis scheme works well for ^^Ti, only in this case, in contrast with the 
previous results for ^^Mg where SU (3) was found to be dominant and SSM 
recessive, in the lower p/-shell SSM is dominant and SU (3) recessive. 

3 Boson-based Applications of the Theory 

In this section, we consider the transition from the U{5) vibrational to the 
SU{3) rotational limit of the interacting boson model using a schematic 
Hamiltonian. The behavior of low-lying energy levels and E2 transition rates 
are studied in detail. The analysis shows that as one moves from the SU (3) 
to the U (5) limit the system transforms from one of maximum deformation 
to one that is maximally triaxial and soft with the onset of strong triaxiality 
occuring around the X(5) critical point. 

Understanding shape phase transitions of a finite many-body system is 
paramount to understanding the system's underlying dynamics. The three 
possible phases that can occur in the interacting boson model for nuclei have 
been classified as U(5), SU(3), and 0(6) [9]. The f/(5) ^ SU{3) transitional 
description of the rare-earth nuclei Nd, Sm, Gd, and Dy was first reported 
in [10], including detailed results for most quantities of physical interest. 
Evidence for coexisting phases at low energy in the spherical-deformed tran- 
sitional nucleus ^^^Sm was also analyzed using C/(5) <-> SU{3) transitional 
theory and the results show that the two phases coexist in a very small re- 
gion of parameter space around the critical point [11]. Recently, since the 
discovery of the X(5) symmetry in this region [12], the spherical to axially 
deformed shape phase transition has attracted further attention [13] . 

In order to take a close look at the U (5) <-> SU (3) shape phase transition. 
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we study transitional patterns of many physical quantities, such as low-lying 
energy levels, isomer shifts, E2 transition rates, and some related quantities 
across the U (5) <-> SU (3) leg of the Casten triangle. No attempt is made to 
relate the results to realistic nuclei; rather, our purpose is to gain a better 
understanding the nature of the U (5) ^ SU (3) transition. In the study, the 
schematic Hamiltonian 



which has been suggested as a suitable form for describing nuclei in this 

region is used, where the parameter c>0, 0<a;<lis the phase parameter, 
f{N) is a linear function of the total number of bosons A^, fis = s^s is the 
number of s bosons, and Q = (s'^d+Ss) — ^(d)d)'^'^\ We note that f{N) = 1 
was used in [11], while f{N) = iN was adopted in [12] and [14]. Also, the 
critical point Xc will be quite different for different choices of the function 
f{N). The Hamiltonian (1) is, up to a constant, equivalent to the one used 
in [10, 11, 12] and [14] with the relation ( = 1 — x. 



In order to diagonalize the Hamiltonian (1), we expand eigenstates of (1) 
in terms of the U{6) D SU{3) D S0{3) basis vectors \N{Xi^)KL) as 



where the quantum number ^ indicating the ,^-th level with angular momen- 
tum quantum number L is introduced, and C^x^^^k expansion coeffi- 
cients. Since the total number of bosons N is fixed for a given nucleus, the 
above eigenstates will also be denoted as \L^;x) in the following, in which 
value of the phase parameter x is explicitly shown. In our calculation, the 
orthonormalization process with respect to the band label K and the phase 
convention for the S'f/(3) D 5'0(3) basis vectors proposed in [15, 16] is used. 
By using analytic expressions for U (6) D SU (3) reduced matrix element of 
s-boson creation or annihilation operator [15] and an algorithm [16, 17] for 
generating the SU{3) D SO (3) Wigner coefficients, the eigen-equation that 
simultaneously determines the eigencncrgy and the corresponding set of the 
expansion coefficients {C^^^^^j^} can be established. 

To explore the transitional patterns, we fix the total number of bosons at 
= 10 and allow the phase parameter x to vary in the closed interval [0, 1]. 




(1) 



(2) 
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Figure 7: Some low-lying energy levels across the transitional region. 
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The functional form of f{N) is chosen to be the same as that used in [12] 
with f{N) = AN unless otherwise specifically noted. Some low-lying energy 
levels as a function of x are shown in Figure 7 from which one can see that 
there is a minimum in the excitation energy around x ~ 0.41 — 0.46 , which 
corresponds to the spherical-deformed shape coexistence region which is also 
referred to as the critical (phase transition) region. It can also be seen that 
the minimum is not exactly the same for all the levels, but all fall within the 
spherical-deformed shape coexistence region. 




Figure 8: Energy ratios i?02 and i?42 as functions of the transitional parameter 

X. 

Two energy ratios, i?02 = Eq+ / E2+ and R42 — / E2+ , are shown as a 
function of x in Figure 8. The ratio i?o2 drops rather precipitously from the 
axially deformed limit, i?02 = 25.338, to the spherical limit, -R02 ~ 2, over 
the range < a; < ~ 0.5 and then remains at the spherical limit value 
for X > 0.5. While the ratio drops rather smoothly as a function of x 
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from the deformed limit, R42 = 10/3 when < a; < ~ 0.3, to the spherical 
limit, i?42 = 2 when ~ 0.6 < x < 1. The sharpest change occurs around the 
critical point Xc ~ 0.46 when the absolute value of the derivative of R42 with 
respect to x reaches the maximal value. 

To show how the transition occurs in the ground state, the amplitudes 
|C(A^)P arc plotted as functions of ,t in Figure 9, which indicate that the most 
rapid changes in these amplitudes also occur within the coexistence region. 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 

X X 



Figure 9: Amplitudes |C(A/i)P of the ground state as functions of the transi- 
tional parameter x. 

In order to explore the exact nature of the critical point with different 
choices for the linear function f{N), overlaps | (0^; x|Og ; a;o) | with = or 
1 as suggested in [12] with f{N) = 0.5N, 2N, AN, and 8N, respectively, 
were calculated. The results are shown in Figure 10, from which it can be 
clearly seen that the critical (crossing) point changes with different choices of 
the function f{N). The larger the f{N) value, the smaller the critical point 
value Xc- This conclusion confirms the early result shown in [11], in which the 
critical point value is very large with Xc = 0.974 corresponding to f{N) = 1, 
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Figure 10: Overlaps of the ground state wavefunction, where the full line 
shows the overlap |(Og;a; = 0|Og;a;)|, and the dotted line shows the overlap 
|(0+;,T = l|0+;a:)|. (a) f{N) = Q.5N; (b) f{N) = 2N- (c) /(TV) = AN- (d) 
f{N) = 8N. 
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while for x, ~ 0.46 when f{N) = AN with = 10 is used [12, 14]. It 
should also be pointed out that the critical point may differ from one excited 
state to the next for the same choice of f{N). Our calculation shows that the 
critical points for Og , 2f, and if are almost the same, but they are somewhat 
different for 2f and 3^. It should also be noted that curves of the overlaps 
\{Lf;x\Lf;xQ)\ with Xq = or 1 even become irregular for 0^, 2^, 4^, and 
O3 . Typical examples of these curves for higher excited states are shown in 
Figure 11, which indicate that the critical point for higher excited states may 
be quite different from that for ground state. This behavior will be seen in 
other observables such as E2 transition rates since the critical behavior of 
the initial and final states may quite different from that of ground state. 




Figure 11: Typical overlaps of two excited states, where the full line shows the 
overlap \{L^;x = 0\L^;x)\, and the dotted line shows the overlap \{L'^;x = 
l\Lp,x)\. 

Various B(E2) values and ratios among the low- lying levels were studied 
using the E2 operator is T{E2) = q2Q, where q2 is the effective charge. 
Some B(E2) values and ratios were found to be sensitive to the shape phase 
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Figure 12: Some B(E2)/g| values as functions of the transitional parameter 
X, where curve a represents B(E2, 23" — > 0+)/g'2, curve b represents B(E2, 
2j 2f)/q2, and curve c represents B(E2, 2j Os")/?!- 
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transition. Figure 12 provides three B(E2) values, B(E2, 2;[ 0+), B(E2, 
23" — i> 2f ), and B(E2, 23" — > O2" ). There is a small peak around x ~ 0.44 for 
B(E2, 2j 0+) and B(E2, 2j ^ 2^), while there is a saddle point around 
X ~ 0.38 for B(E2, 23 — O2 )• As seen previously, the peak and saddle points 
in these B(E2) values are different from the critical point of the ground state 
due to different critical behavior of the excited states. Figure 13 provides six 
B(E2) ratios, B(E2; 2^ ^ 0+)/ B(E2; 2^ ^ 2^), B(E2; 2^ ^ 2^)/ B(E2; 
2+ ^ 0+), B(E2; 2+ ^ 0+)/ B(E2; 2+ ^ 0+), B(E2; 2+ ^ 0+)/ B(E2; 
2+ 0+), B(E2; 2+ ^ 2+)/ B(E2; 2+ ^ 0+), and B(E2; 3+ ^ 2+)/ B(E2; 
3^ Af ). These ratios all undergo noticeable changes within the coexistence 
region. The most distinctive signature is shown in B(E2; 23" — >■ Of )/ B(E2; 
23 — >• O2 ), in which there is a giant peak around x ~ 0.38. 

In order to study shape of nuclei around the critical point, we use the 
relation between the Bohr variables {/3, 7) of the collective model and the 
(A,/i) labels that define the irreducible representation of the SU{3) [18]. In 
this algebraic approach, the Bohr variables can be expressed as a functional 
of SU{3) Casimir operator of the second order P = Pq\JC2{SU{3)) + 3, and 
7 = tan^^ ( ^^^^^} ) , where Pq can be taken as a constant and A and fi should 
be regarded as operators, of which the results are usual A and /j, values when 
they are applied onto the basis vector of SU (3). 

Using these definitions for (3 and 7, we can calculate expectation values, 
P — (Og ; x|;9|0g ; x) and 7 = (Og ; x|7|0g ; x) in the ground state, and the 
corresponding root mean square deviations A(/3) and A(7). These quantities 
can be used to display shape uncertainty as a function of x. It is obvious 
that A(/3) and A (7) are zero when nucleus is spherically deformed with x — 
indicating there is a definite shape, while the shape becomes less well defined 
as X moves away from zero. The values of A(/9) and A(7) are not small 
in the spherical limit (a; = 1) due to the quadrupole vibration being not 
negligible. It can be seen from Figure 14 that there are also obvious changes 
in P and 7 in the critical region around x ~ 0.46. Over the whole range, the 
magnitude of 7 change is small. In contrast, A(/3) reaches a maximum value 
around x ~ 0.5 which deviates a little from the critical point 0.46 of the 
ground state, but is still near the critical region. This distinctive signature 
shows that nucleus is the softest in the critical region. This fact can help us 
to understand why there is a saddle region in most excited levels. Since a 
nucleus in this region is comparatively soft, its shape can be changed easily 
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Figure 13: Some B(E2) ratios as functions of the transitional parameter 
where curve a represents B(E2; 2^ — > 0+)/ B(E2; 23 — > 2^), curve b 
represents B(E2; 2^ 2|)/ B(E2; 2f 0+), curve c represents B(E2; 
'^t 0^)/ B(E2; 2t 0+), curve d represents B(E2; 2^ ^ 0^)/ B(E2; 
23 ^ 0^), curve e represents B(E2; 2^ 2^)/ B(E2; 2^ 0+), and curve 
f represents B(E2; 3^^ ^ 2t)/ B(E2; 3^^ ^ 4^"). 
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with very little energy. Hence, nuclei in this region can easily be excited, 
which results in relatively smaller energy gaps in this soft critical region. 




Figure 14: The ground state expectation values of Bohr variables P and 7 
and the corresponding root mean square deviations A(/9) and A (7), where 
the full line indicates the expectation value f3/f3Q or 7, while the dotted lines 
show the corresponding root mean square deviations ±A(/3) or ±A(7). 

In summary, transitional patterns from the vibrational, U{5), to the rota- 
tional, SU{3), limit of the interacting boson model with a schematic Hamil- 
tonian have been studied. The transitional behavior of low-lying energy lev- 
els, isomer shifts, E2 transition rates, and related quantities over the whole 
U (5) <-> SU (3) transitional region were explored. The results show that there 
are many distinctive signatures in the energy levels, wavefunctions, isomer 
shift, BE2 values and ratios, and expectation value of Bohr variables near 
the critical point. Generally speaking, critical the behavior of excited states 
are different from that of ground state, which may lead to different critical 
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point for some physical quantities that hnk the two. In comparison with the 
result shown in [20], a nucleus with X(5) symmetry can approximately be 
described by the U (5) <-> SU (3) transitional theory within the critical region. 
Our analysis also shows that shapes of nucleus in the critical region is not 
well-dfined; that is, a nucleus with X(5) symmetry is soft. 

4 Conclusions 

Future research may provide justification for an extension of the theory to 
multi-mode shell-model calculations. For example, an immediate extension 
of the current scheme could use eigenvectors of the pairing interaction [19] 
within an S'p(4) algebraic approach to the nuclear structure [21] in addition 
to collective SU (3) configurations and spherical shell model states. Alterna- 
tively, Hamiltonian driven basis sets could be considered. For example, the 
method could use eigenstates of near closed shell nuclei obtained from a full 
shell-model calculation to form Hamiltonian driven J-pair states for mid-shell 
nuclei [22]. This would mimic the Interacting Boson Model (IBM) [23] and 
the so-called broken-pair theory [22]. Likewise, the three exact limits of the 
IBM [24] can be considered to comprise a three-mode system. Nonetheless, 
the real benefit of the mixed-mode approach is expected when the spaces 
encountered are too large to allow for exact calculations. 

This work was supported by the U.S. National Science Foundation (0140300) 
and by the Natural Science Foundation of China (10175031). 
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